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We describe the transverse beam distribution in particle accelerators within the controlled, 
stochastic dynamical scheme of the Stochastic Mechanics (SM) which produces time reversal in- 
variant diffusion processes. This leads to a linearized theory summarized in a Shchrodinger-like 
(S-^) equation. The space charge effects have been introduced in a recent paper l] by coupling this 
S-£ equation with the Maxwell equations. We analyze the space charge effects to understand how 
the dynamics produces the actual beam distributions, and in particular we show how the stationary, 
self-consistent solutions axe related to the (external, and space-charge) potentials both when we 
suppose that the external field is harmonic {constant focusing), and when we a priori prescribe the 
shape of the stationary solution. We then proceed to discuss a few new ideas by introducing the 
generalized Student distributions, namely non-Gaussian, Levy infinitely divisible (but not stable) 
distributions. We will discuss this idea from two different standpoints: (a) first by supposing that 
the stationary distribution of our (Wiener powered) SM model is a Student distribution; (b) by 
supposing that our model is based on a (non-Gaussian) Levy process whose increments are Student 
distributed. We show that in the case (a) the longer tails of the power decay of the Student laws, 
and in the case (b) the discontinuities of the Levy-Student process can well account for the rare 
escape of particles from the beam core, and hence for the formation of a halo in intense beams. 

PACS numbers; 02.50.Ey, 05.40.Fb, 29.27.Bd, 41.75.Lx 



I. INTRODUCTION 

In high intensity beams of charged particles, proposed 
in recent years for a wide variety of accelerator-related 
applications, it is very important to keep at low level the 
beam loss to the wall of the beam pipe, since even small 
fractional losses in a high-current machine can cause ex- 
ceedingly high levels of radioactivation. It is now widely 
believed that one of the relevant mechanisms for these 
losses is the formation of a low intensity beam halo more 
or less far from the core. These halos have been observed 
or studied in experiments Q , and have also been sub- 
jected to an extensive simulation analysis 0. For the 
next generation of high intensity machines it is however 
still necessary to obtain a more quantitative understand- 
ing not only of the physics of the halo, but also of the 
beam transverse distribution in general jOj. In fact "be- 
cause there is not a consensus about its definition, halo 
remains an imprecise term" Q so that several proposals 
have been put forward for its description. 

The charged particle beams are usually described in 
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terms of classical, deterministic dynamical systems. The 
standard model is that of a coUisionless plasma where the 
corresponding dynamics is embodied in a suitable phase 
space (see for example Q). In this framework the beam 
is studied by means of the particle-in-core {pic) model 
and the simulations show that the instabilities due to a 
parametric resonance can allow the particles to escape 
from the core with consequent halo formation |5| |^ |3 ■ 
The present paper takes a different approach: it follows 
the idea that the particle trajectories are samples of a 
stochastic process, rather than usual deterministic (dif- 
ferentiable) trajectories. In the usual dynamical mod- 
els there is a particle probability distribution obeying 
the Vlasov equation, and its evolution is Liouvillian in 
the sense that the origin of the randomness is just in 
the initial conditions: along the time evolution, which is 
supposed to be deterministic, there is no new source of 
uncertainty. It is the non linear character of the equa- 
tions which produces the possible unpredictable charac- 
ter of the trajectories. On the other hand in our model 
the trajectories are replaced by stochastic processes since 
the time evolution is supposed to be randomly perturbed 
even after the initial time. It is open to discussion which 
one of these two description is more realistic; in particular 
we should ask if the mutual interactions among the beam 
particles look like random collisions, or rather like con- 
tinuous deterministic interactions. In the opinion of the 
authors, however, a plasma (with collisions) described in 
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terms of controlled stochastic processes is a good candi- 
date to explain the rare escape of particles from a quasi- 
stable beam core by statistically taking into account the 
random inter-particle interactions that can not be de- 
scribed in detail. Of course the idea of a stochastic ap- 
proach is hardly new , but there are several different 
ways to implement it. 

First of all let us remark that the system we want to 
describe is endowed with some measure of invariance un- 
der time reversal, since the external fields act to keep it 
in a quasi-stationary non diffusive state despite the re- 
pulsive electro-magnetic (e.m.) interactions among the 
constituent particles. However, a widespread misconcep- 
tion notwithstanding, a theory of stochastic processes not 
always describe irreversible systems: the addition of a 
dynamics to a stochastic kinematics can in fact ascribe 
a measure of time reversal invariance also to a stochas- 
tic system The standard way to build a stochastic 
dynamical system is to modify the phase space dynamics 
by adding a Wiener noise B(t) to the momentum equa- 
tion only, so that the usual relations between position 
and velocity is preserved: 

mdQ(t) = P{t)dt, 

dP{t) ^ F{t)dt + pdB{t) . 

In this way we get a derivable, but not Markovian po- 
sition process Q(i). The standard example of this ap- 
proach is that of a Brownian motion in a force field de- 
scribed by an Ornstein-Uhlenbeck system of stochastic 
differential equations (SDE) ^Jl- Alternatively we can 
add a Wiener noise W(t) with diffusion coefficient D to 
the position equation: 

dQit) = v^+^{Q{t),t)dt + ^/DdW{t). 

and get a Markovian, but not derivable Q(t). In this way 
the stochastic system is also reduced to a single SDE 
since we are obliged to drop the second (momentum) 
equation: in fact now Q(t) is no more derivable. The 
standard example of this reduction is the Smoluchowski 
approximation of the Ornstein-Uhlenbeck process in the 
overdamped case pTl |. As a consequence we will work 
only in a configuration, and not in a phase space; but 
this does not prevent us from introducing a dynamics - 
as we will show in the Secti on l ll Al - either by generalizing 
the Newton equations Illl.[l3l , or by means of a stochas- 
tic variational principle |l3| | . Remark that in this scheme 
the forward velocity V(_|_) (r, t) can no more be an a priori 
given field: rather it now plays the role of a new dy- 
namical variable of our system. This second scheme, the 
Stochastic mechanics (SM), is universally known for its 
original application to the problem of building a classical 
stochastic model for Quantum Mechanics (QM), but in 
fact it is a very general model which is suitable for a large 
number of stochastic dynamical systems |0, We 
will also see in the Section FlI Al that from the stochastic 
variational principles two coupled equations are derived 
which are equivalent to a Schrodinger-like (S-€) differen- 
tial equation: in this sense we will speak of quantum-like 



(Q-i) systems, in analogy with other recent researches 
on this subject [H In fact the SM can be used 

to describe every stochastic dynamical system satis fying 
fairly general conditions: it is known since longtime |l7l| . 
for example, that for any given diffusion there is a cor- 
respondence between diffusion processes and solutions of 
S-^ equations where the Hamiltonians come in general 
from suitable vector potentials. Under some regularity 
conditions this correspondence is seen to be one-to-one. 
The usual Schrodinger equation, and hence QM, is recov- 
ered when the diffusion coefficient coincides with S/2m, 
namely is connected to the Planck constant. However we 
are interested here not in a stochastic model of QM, but 
in the description of particle beams. 

In the present paper we intend to widen the scope 
of our SM model by introducing the idea that an im- 
portant role for the beam dynamics can be played by 
non-Gaussian Levy distributions. In fact these distribu- 
tions enjoyed a widespread popularity in the recent years 
because of their multifaceted possible applications to a 
large set of problems from the statistical mechanics to 
the mathematical finance (see for example 0, 0| and 
references quoted therein). In particular the so called 
stable laws (see Section ^ are used in a large number of 
instances, as for example in the definition of the so-called 
Levy flights. Our research is instead focused on a family 
of non-Gaussian Levy laws which are infinitely divisible 
but not stable: the generalize Student laws. As will be 
discussed later this will allow us to overcome - without 
resorting to the trick of the truncated laws - the prob- 
lems raised by the fact that the stable non-Gaussian laws 
always have divergent variances: a feature which is not 
realistic to ascribe to most real systems. It is possible 
to show indeed that by suitably choosing the parame- 
ters of the Student laws we can have distributions with 
finite variance, and approximating the Gaussian law as 
well as we want. On the other hand the infinitely divisi- 
ble character of these laws is all that is required to build 
a stationary, stochastically continuous Markov process 
with independent increments, namely the Levy process 
that we propose to use to represent the evolution of our 
particle beam. 

Of course it is not always mathematically easy to deal 
with the infinitely divisible processes, but we will show 
that at least in two respects they will help us to have some 
further insight in the beam dynamics. First of all we use 
the Student distributions in the framework of the tradi- 
tional SM where the randomness of process is supplied by 
a Gaussian Wiener noise: here we examine the features 
of the self-consistent potentials which can produce a Stu- 
dent distribution as stationary transverse distribution of 
a particle beam. In this instance the focus of our research 
is on the increase of the probability of finding the par- 
ticles at a great distance from the beam core. Then we 
pass to the definition of a true Levy-Student process, and 
we show with a few simulations that these processes can 
help to explain how a particle can be expelled from the 
bunch because of some kind of hard collision. In fact the 
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trajectories of our Levy-Student process show the typi- 
cal jumps of the non-Gaussian Levy processes: a feature 
that we propose to use as a model for the halo formation. 
It is worth remarking that, albeit the more recent empir- 
ical data about halos still not accurate enough to 
distinguish between the suggested distributions and the 
usual Gaussian ones, our conjecture on the role of Stu- 
dent laws in the transverse beam dynamics has recently 
found a first confirmation poj | in numerical simulations 
showing how these laws are well suited to describe the 
statistics of the random features of the particle paths. 

In a few previous papers we connected the (trans- 
verse) r.m.s. emittance to the characteristic microscopic 
scale and to the total number of the particles in a bunch, 
and implemented a few techniques of active control for 
the dynamics of the beam. In this paper we first of all 
review the theoretical basis [3, |^ of the proposed model: 
in the Section |n] we define our SM model with emphasis 
added on the potentials which control the beam dynam- 
ics and on the possible non stationary solutions of this 
model |25|. In the Section lllll we review our analysis of 
the self-consistent, space charge effects due to the e.m. 
interaction among the particles, adding a few new results 
and comments. In the Section IIVI we then discuss the 
idea that the laws ruling the transverse distribution 
of particle beams are non-Gaussian, infinitely divisible. 
Levy laws as the generalized Student laws. In particular 
we analyze the behavior of our usual SM model under 
the hypothesis that the stationary transverse distribu- 
tion is a Student law. Finally in the Section we study 
the possibility of extending our SM model to Levy pro- 
cesses whose increments are distributed according to the 
Student law. We think in particular that the presence 
of isolated jumps in the trajectories can help to build a 
realistic model for the possible formation of halos in the 
particle beams. We end the paper with a few conclusive 
remarks. 



II. STOCHASTIC BEAM DYNAMICS 

A. Stochastic mechanics 

First of all we introduce the stochastic process per- 
formed by a representative particle that oscillates around 
the closed ideal orbit in a particle accelerator. We con- 
sider the 3-dimensional (3-DIM) diffusion process Q(i), 
taking the values r, which describes the position of the 
representative particle and whose probability density is 
proportional to the particle density of the bunch. As 
stated in the Section ^ the evolution of this process is 
ruled by the Ito stochastic differential equation (SDE) 

dCl{t)^v^+■,{Q{t),t)dt + ^/DdW{t), (1) 

where V(+)(r, i) is the forward velocity, and dW{t) = 
W(t -f dt) — W(i) is the increment process of a standard 
Wiener noise W(i); as it is well known this increment 



process is gaussian with law Af{0,ldt), where I is the 
3x3 identity matrix. Finally the diffusion coefficient 
D is supposed to be constant: the quantity a = 2mD, 
which has the dimensions of an action, will be later con- 
nected to the characteristic transverse emittance of the 
beam. The equation defines the random kinematics 
performed by the particle, and replaces the usual deter- 
ministic kinematics 

dq(t) =v(q(t),t)dt (2) 

where q(i) is just the trajectory in the 3-DIM space. 

To counteract the dissipation due to this stochastic 
kinematics, a dynamics must be independently added. 
In SM we do not have a phase space: our description 
is entirely in a 3-DIM configuration space. This means 
in particular that the dynamics is not introduced in a 
Hamiltonian way, but by means of a suitable stochastic 
least action principle obtained as a generalization of 
the variational principle of classical mechanics. In the 
following we will briefly review the main results, refer- 
ring for details to the references [Til, [ill IT3j | . Given the 
SDE JQ), we consider the probability density function 
(pdf) p{r,t) associated to the diffusion Q(t) so that, be- 
sides the forward velocity V(+) (r, t), we can now define a 
backward velocity 

V(_)(r,0=V(+)(r,t)-2i^^M. (3) 

We can then introduce also the current and the osmotic 
velocity fields, defined as: 

v="W+"'-^ u="W(-> =D^. (4) 
2 2 p 

Here v represents the velocity field of the density, while 
u is of intrinsic stochastic nature and is a measure of the 
non differentiability of the stochastic trajectories. 

A first consequence of the stochastic generalization of 
the least action principle 1 1 , Ts] is that the current ve- 
locity takes the following irrotational form: 

TOv(r,t) = VS'(r,i) , (5) 

while the Lagrange equations of motion for the density p 
and for the current velocity v are the continuity equation 
associated to every stochastic process 

dtp = -V • (pv) , (6) 

and a dynamical equation 

d,S+^w'^2mD^^^ + V{r,t)^0, (7) 
2 ^/P 

which characterizes our particular class of time-reversal 
invariant diffusions (Nelson processes). The last equa- 
tion has the same form of the Hamilton-Jacobi- 
Madelung (HJM) equation, originally introduced in the 
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hydrodynamic description of quantum mechanics by 
Madelung 23]. Since (|SJ) holds, the two equations © 
and can be put in the foUowing form 

dtp = --V-(pV5) (8) 
m 

dtS = -^VS^ + 2mD^^^ -V{r,t) (9) 
2m 

which now constitutes a coupled, non linear system of 
partial differential equations for the pair (p, S) which 
completely determines the state of our beam. On the 
other hand, because of (|SJ), this state is equivalently given 
by the pair(/9, v). 

It can also be shown by simple substitution from Q 
that © is equivalent to the standard Fokker-Planck 
(FP) equation 

5tp = -V-[v(+)p] + i?VV (10) 

formally associated to the Ito equation In fact also 
the HJM equation Q can be cast in a form based on 
V(^) rather than on v, namely 

dtS = -yvf+)+mDv(+)Vln/ 

+mD^V^liif -V (11) 

where / is a dimensionless density defined by 

p(r,t) = C/(r,t) (12) 

where C is a dimensional constant. On the other hand, 
from Q and Q , we know that also the forward velocity 
V(+) is irrotational: 

V(+)(r,t)=VM/(r,t), (13) 

and that by taking (O into account the functions W and 
S are connected by the relation 

S{r, t) = mW{v, t) - mD In /(r, t) - 6{t) (14) 

where 9 is an arbitrary function of t only. 

The time-reversal invariance is now made possible ^3 
by the fact that the forward drift velocity V(_|.) (r, t) is no 
more an a priori given field, as is usual for the diffusion 
processes of the Langevin type; instead it is dynamically 
determined at any instant of time, starting by an initial 
condition, through the HJM evolution equation (Tj) . It is 
finally important to remark that, introducing the repre- 
sentation [23 

*(r,i) = v/;^e'«('-'*)/", (15) 

(with a — 2mD) the coupled equations ^ and ® are 
made equivalent to a single linear equation of the form 
of the Schrodinger equation, with the Planck action con- 
stant replaced by a: 

2 

iadt"^ ^ -^V^i; + V'^ . (16) 
2m 



We will refer to it as a Schrodinger-like (S-£) equa- 
tion: clearly (|16|l has not the same meaning as the usual 
Schrodinger equation; this would be true only if a — h, 
while in general a is not an universal constant, and it 
is rather a quantity characteristic of the system under 
consideration (in our case the particle beam). In fact a 
turns out to be of the order of magnitude of the beam 
emittance, a quantity which - in formal analogy with h 
- has the dimensions of an action and gives a measure of 
the position/momentum uncertainty product for the sys- 
tem. Thus the SM model of our beam, as incorporated 
in the phenomcnological Schrodinger equation ]l<i\i . while 
keeping a few features reminiscent of the QM, is in fact 
a deeply different theory. 



B. Controlled distributions 

We have introduced the equations that in the SM 
model are supposed to describe the dynamical behavior 
of the beam: we now briefly sum up a general procedure, 
already exploited in previous papers 2L .24] , to control 
the dynamics of our systems. Let us suppose that the 
pdf p(r, t) be given all along its time evolution: think 
in particular either to a stationary state, or to an en- 
gineered evolution from some initial pdf toward a final 
state with suitable characteristics. We know that the 
FP equation 1)10(1 must be satisfied, for the given p, by 
some forward velocity field V(+) (r, t). Since also the equa- 
tion l|13() must hold, we are first of all required to find an 
irrotational V(+) which satisfies the FP equation ((10|l for 
the given p. Wc then take into account also the dynami- 
cal equation (|11|) : since p and v^^) (and hence / and W) 
are now fixed and satisfy ((l()|l . the equation ((11(1 plays 
the role of a constraint defining a controlling potential V 
when we also take into account the equation ((14(1 . We 
list here the potentials associated to the three particular 
cases analyzed in the previous papers. 

In the 1-DIM case with given dimensionless pdf f{x, t) 
and a < X <b (a and b can be infinite) we easily get 

V{x, t) = mD'^ dl In / + mD {dt In / -t- v^+-)d^ In /) 

~—vl^^~mj dtV(+){x' ,t)dx' + e (18) 

For a 3-DIM system with cylindrical symmetry around 
the z-axis (the beam axis) , if we denote with (r, z) 
the cylindrical coordinates, and if we suppose that p{r, t) 
depends only on r and t, and that V(_|_) = (r, t) f is 
radially directed with modulus depending only on r and 
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i, we have 

^(+)(^^)-^%^-r77^ / dtp{r',tydr' (19) 



p(r, t) rp{r, t) Jq 



mD 



V{r, t) = dr{rdr In /) + mD {dt In / + V(^+^dr In /) 



- 2" «(+) - "1 dtvi^+) (j' ,t)dr' + 6 (20) 

Finally in the 3-DIM stationary case the pdf p{r) is in- 
dependent from t. This greatly simplifies our formulas 
and, by requiring that 9{t) ~ E he constant, namely 
that 9{t) = Et, we get 



v(+)(r) = D 



Vp(r) 



p(r) 

V{r) = E + 2mD 



2 v^VP 



VP 



(21) 
(22) 



Of course in this context the constant E will be chosen 
by fixing the zero of the potential energy. Let us remark 
finally that in this stationary case the phenomenological 
wave function (|15|l takes the form 



*(r,<) = 
typical of the stationary states. 



'iEtja 



C. Non stationary distributions 

In the following we will be mainly concerned with sta- 
tionary distributions, but in a few previous paper we 
treated also non stationary problems. For instance, if 
we consider the stationary, ground state pdf (without 
nodes) po(r) of a suitable potential, and if we calculate 
V(^)(r) and write down the corresponding FP equation, 
it is possible to show (see the general proof in a few pre- 
vious papers ISEHEBl) that, /Oo(r) will play the role of 
an attractor for every other distribution (non extremal 
with respect to a stochastic minimal action principle) . If 
the accelerator beam is ruled by such an equation, this 
would imply that the halo can not simply be wiped out 
by scraping away the particles that come out of the bunch 
core: in fact they simply will keep going out in the halo 
until the equilibrium is reached again since the distribu- 
tion po(r) is a stable attractor. 

In a recent paper [23 we gave an estimate of the time 
required for the relaxation of non extremal pdf 's toward 
the equilibrium distribution. This is an interesting test 
for our model since this relaxation time is fixed once the 
form of the forward velocity field is given; this is in turn 
fixed when the form of the halo distribution is given as 
in the reference ^J, and one could check if the estimate 
is in agreement with possible observed times. In par- 
ticular we estimated that in typical conditions all the 
non-stationary solutions of this FP equation will be at- 
tracted toward po with a relaxation time of the order of 
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A different non stationary problem also discussed in 
previous papers |^ consists in the analysis of some 
particular time evolution of the process with the aim of 
finding the dynamics that control it. For instance we 
studied the possible evolutions which start from a pdf 
with halo and evolve toward a halo-free pdf: this would 
allow us to find the dynamics that we are requested to 
apply in order to achieve this result. If for simplicity the 
overall process is supposed to be an Ornstein-Uhlenbeck 
process, the transition pdf would be completely known 
and all the result can be exactly calculated through the 
Chapman-Kolmogorov equation by supposing suitable 
shapes for the initial and final distributions. Then a di- 
rect application of H18|l allows us to calculate the control 
potential corresponding to this evolution. For the sake 
of brevity we do not give the analytical form of this po- 
tential and refer to the quoted papers for further details. 



III. SELF-CONSISTENT EQUATIONS 

A. Space charge interaction 

In QM a system of N particles is described by a wave 
function in a 3iV-DIM configuration space. On the other 
hand in our SM scheme a normalized |5'(r,t)p, func- 
tion of only three space coordinates r = {x,y,z}, plays 
the role of the pdf of a Nelson process. In a first ap- 
proximation we will consider this iV-particle system as 
a pure ensemble: as a consequence we will not introduce 
a SA^-DIM configuration space, since N\^{Y^t)\'^ (Pr in 
the 3-DIM space will play the role of the number of par- 
ticles in a small neighborhood of r. However, since our 
system of N charged particles is not a pure ensemble 
due to their mutual e.m. interaction, in a further mean 
field approximation we will take into account the so called 
space charge effects: more precisely we will couple our S- 
l equation with the Maxwell equations describing both 
the external and the space charge e.m. fields, and we will 
get in the end a non linear system of coupled differential 
equations. 

In our model a single, charged particle embedded in a 
beam and feeling both an external, and a space charge 
potential is first of all described by a S-£ equation 

ia9f^'(r,i) = ^*(r,t) , 

where vE'(r, t) is our wave function, a a coefficient with the 
dimensions of an action which is a constant depending on 
the beam characteristics, and iJ is a suitable Hamiltonian 
operator. If ^ is properly normalized and if N is the 
number of particles with individual charge goi the space 
charge density and the electrical current density are 

Psc{v,t) = iVqo|*(r,OI', (23) 
j,,(r,i) = iV<7o-3{**(r,i)V*(r,i)}. (24) 

TO 

Hence our particles in the beam will feel both an electri- 
cal and a magnetic interaction and we will be obliged to 
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couple the S-£ equation with the equations of the vector 
and scalar potentials associated to this electro-magnetic 
field. 

The e.m. potentials (A^c, ^sc) of the space charge fields 
obeying the gauge condition 



V • A,,(r,t) + —dt^scir,t) = 0, 



(25) 



must satisfy the wave equations 
V2<i>,,(r,t)-la2^,,(r,t) 



-^ojsc(r,i) (26) 
P.c(r,<) ^27) 



On the other hand, for our particle in the beam the 
e.m. field is the superposition of the space charge poten- 
tial (Asc, ^sc), and of the external potentials (Ae,$e)- 
Hence (see for example j23|, chapter XV) our S-£ equa- 
tion takes the form 
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ioSl (A^c + Ae) 



-'Zo(*.c + $e)«' (28) 



2m 



It is apparent now that (|2SJ|, and (|^ con- 

stitute a self-consistent system of non linear differen- 
tial equations for the fields 5*, A^c and $sc coupled 
through and it^ . 

If we then consider stationary wave functions 



(29) 



where E is the energy of the particle, and take Ae = for 
the external interaction, passing to the potential energies 

Ve(r) = go*e(r) , Vsc[r:) = go*sc(r) , 

our system is reduced to only two coupled, non linear 
equations for the pair (-0, t^sc), namely 



2m 



(K + Vsc - E)^ , 



(30) 
(31) 



B. Cylindrical symmetry 



We suppose now that the longitudinal motion along 
the z~axis is both decoupled from the transverse motion 
in the x, y-plane, and free with constant momentum pz, 
and velocity bz = bo ^ bx,by. Moreover we suppose that 
the beam particles will be confined in a cylindrical packet 
of length L, so that by the imposing periodic boundary 
conditions we will quantize the longitudinal momentum 



2A;7rQ; 
L 



0,±1,±2,.. 



As a consequence our wave functions will take the form 



and our equations H30() and H31|l become 



{dl + dl)Vs. 



Nqf 
Le 



(0 



(32) 

(33) 
(34) 



where M = N / L is the number of particles per unit 
length, and Et ~ E — p1/2m is the energy of the trans- 
verse motion. If finally our system has a cylindrical sym- 
metry around the z axis, namely if - in the cylindrical 
coordinate system {r, ip, z} (r^ = x'^ + y^) - our poten- 
tials depend only on r, then we can separate the variables 
with x(x,?/) = u{r)<^{ip), the angular eigenfunctions are 



7^ 



^ = 0,±1,±2,, 



and for £ — the equations become 

^,2 



2m 



V" 



u 
r 

V 

^ Si 

r 



(14 + Vsc - Et) u 

^ll 2 

— u 

27reo 



(35) 

(36) 
(37) 



with the following radial normalization 

f + C30 

ru (r) dr = \ . 



Remark that now we are reduced to a system of ordinary 
differential equations. 



C. Dimensionless formulation 

To eliminate the physical dimensions one introduces 
two quantities and A which are respectively an energy 
and a length. Then, by means of the dimensionless quan- 
tities 
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A 77 



Afgl 
27reo77 



(perveance) 



w{s) — A u{Xs) 



VscjXs) K(A,s) 

V{S) = , Ve[S) = 



V ' V 

the equations H36|) and H37|l take the form 

sw"{s)+w'{s) = [ve{s) + v{s) ~ P] s w{s) (38) 
sv"{s)+v'{s) = -^sw^{s) (39) 

The usual choice for the dimensional constants is 



r] — mbl , A = 



mb[)^/2 ' 



(40) 
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where Bq is the longitudinal velocity of the beam. We can 
now look at our equations in two different ways. First of 
all we can suppose that Ve is a given external potential: 
in this case our aim is to solve the equations for the two 
unknowns w (radial particle distribution) and v (space 
charge potential energy). However in general no simple 
analytical solution of this problem is at present available 
for the usual forms of the external potential Ve '■ there are 
not even solutions playing the same role played by the 
Kapchinskij-Vladimirskij (KV) distribution in the usual 
models. This phase space distribution - which is simple 
and self-consistent in the usual dynamical models - leads 
to an uniform transverse space distribution of the beam, 
and is a stationary solution of the Vlasov equation with a 
harmonic potential. Moreover its space charge potential 
calculated from the Poisson equation is still harmonic. 
Instead in the SM model the uniform distributions are 
not solutions of the stationary Schrodinger equation, and 
we know no simple stationary distribution connected to 
the harmonic potential as the KV. Even the gaussian dis- 
tributions - later discussed in this paper - can not play 
the same role: they are solutions connected with an exter- 
nal harmonic potential, but their space charge potential 
calculated from the Poisson equation is not harmonic. 

Alternatively we can assume as known a given distribu- 
tion w, and solve our equations to find both the external 
and the space charge self-consistent potential energies Ve 
and V. In this second form the problem is more simple, 
and analytical solutions are available. We adopted the 
first standpoint in a few previous papers Q where we 
numerically solved the equations (|38|l and H39|l : here we 
will rather elaborate a few new ideas about the second 
one. To this end it is important to remark that the space 
charge potential energy 



vis) = 



dy 

y 



xuP'ix) dx 



(41) 



is always a solution of the Poisson equation %V^\ satis- 
fying the conditions w(0"'") = v'iO'^) = 0. On the other 
hand, by substituting (|1T|| in the first equation we 
readily obtain also the self-consistent form of the exter- 
nal potential energy 



Ve{s) 
wo(s) 



Wo(s)-|-f / — / xw^{x)dx. 

Jo y Jo 
w"is) Iw'(s) 



w{s) s w{s) 



-/3 



(42) 
(43) 



where vo(s) is the potential that we would have with- 
out space charge (^ = 0), while the second part in the 
external potential (j42|l exactly compensate for the space 
charge potential. 



D. Constant focusing 

Let us suppose now that the transverse external poten- 
tial Ve{r) is a cylindrically symmetric, harmonic potential 



with a proper frequency uj (constant focusing) , and let we 
also introduce the characteristic length 



2muj 



which will represents a measure of the transverse disper- 
sion of the beam. In cylindrical coordinates {r, ip\ in the 
transverse plane our potential energy is 



Ve{r) 



8ma 



(44) 



so that the corresponding 2-DIM S~£ equation without 
space charge (zero perveance) would have as lowest eigen- 
value Eq — auj, and as ground state wave function 



Xoo{r, if) 



uo(r) 



(45) 



Of course the self-consistent solution would be different 
if there is a space charge (non zero perveance). To find 
this solution one introduces the so called phase advance 

I Lo a 
Xq bo 2m6oc^ 

(Ao is a length) and, with the constants 14U|) . the dimen- 
sionless form of the harmonic potential H44|l 



Veis) 



Veir) 



UJ 



mbl hi 



2A§ 4A2 m2 62 



2 2 2 

s = 7 s 



2Aoto6o 2mbQ 



^0 



As a consequence the equations (|38f) and (|39|l become 

sw"{s)+w'{s) = [-1^ s^ + v{s)- (3\sw{s) (46) 



sv"{s) + v'{s) = -£,sw'^{s) 



(47) 



These equations are now a coupled, non linear system 
which must be numerically solved since we do not know 
simple self-consistent solutions of the form of the KV dis- 
tribution. In reference we extensively analyzed these 
numerical solutions and we refer to this paper for details. 
In fact in jlj] there was a small difference with respect to 
what has been presented here. The form of the equations 
to solve is the same, but the dimensionless formulation 
was achieved by means of two numerical constants differ- 
ent from (|4U|) and drawn from the characteristics of the 
transverse harmonic oscillator force: 

Then the dimensionless quantities have a different numer- 
ical value and the dimensionless equations 136|l and 137|l 
take the form 



sw"{s) + w'{s) 
sv"{s) + v'{s) 



[s^ + v{s) ~ (3]sw{s) (49) 
-S,sw^{s) (50) 
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since now 7 = 1. In any case the equations (|46() and H47() 
can easily be turned into the equations (|49|l and (|50|l . and 
vice versa, by means of simple transformations through 
the parameter 7 which turns out to be at the same time 
the ratio of the energy constants, and that of the squared 
length constants. As a consequence in the following we 
will always use the system (|49|) . (|50f) . with the advantage 
of simply putting 7 = 1 in the model. 
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IV. SELF-CONSISTENT POTENTIALS 



A. Gaussian transverse distributions 



10 



-50 



In the SM model it is possible to numerically integrate 
the Schrodinger-Poisson system 138() and (|39|l with a 
given external potential and calculate the self-consistent 
distributions and their space charge potentials P|. On 
the other hand, if we fix a particular distribution, it is 
always possible to exactly calculate from these equations 
the external and space charge potential giving rise to that 
distribution. When we adopt this second alternative ap- 
proach and we take as given the form of the distribution 
w{s), the unknowns in the equations (|38|l and (|39|l are 
the two potential energies v(s) and Ve{s). In this case 
we only need to calculate the expressions (|41|l and H42|) 
in terms of the given distribution w{s). Of course if we 
take an arbitrary w{s) we will not get any simple and 
meaningful form for the external potential Ve{s); and on 
the other hand to guess the right form of w{s) giving 
rise, for instance, exactly to a harmonic potential (|44|l as 
external potential would be tantamount to solve (|42|l as 
an integro-differential equation for a given external po- 
tential. However in a few explicit cases the results are 
quite simple and interesting. 

Let us take as first example of a stationary wave func- 
tion that of the ground state uo{r) of the harmonic oscil- 
lator with zero perveance given in 1451) . Its dimensionless 
representation is: 

w;(s) = \/2e"^'/^ /3 = 2 [Et = ato) ; (51) 

which is also apparently normalized. We now want to cal- 
culate both the external and the space charge potentials 
that produce (|51|l as stationary wave function for 149|l 
and Ip1|. From and (|^ we then have 



w"{s) 1 w'{s) 



+ /3 = vois) 



w{s) s w{s) 

r w-^{x)xdx = i nog(s2) + C-Ei(-s2)l 
Jo y Jo 2 L 



where 



0.577 is the Euler constant and 



Ei(x) 



t 



dt, 



X <0 



is the exponential-integral function, and hence we imme- 
diately get (see also FIG.^ 



FIG. 1: The dimensionless potentials v{s) (thin line), vo{s) = 
(dashed line) and Ve{s) (thick line). They reproduce re- 
spectively equations I52II . 11531 and 15411 for ^ = 20 (see refer- 
ence Q for this value). When the external potential is Ve{s) 
the self-consistent wave function coincides with that of a sim- 
ple harmonic oscillator for zero perveance 15H . 



v{s) 

Ve{s) 



-i[log(s^) + C-Ei(-52)] 



[log(s2) + C - Ei(-s2)] 



(52) 
(53) 
(54) 



In a sense the meaning of the equations H41(l . I|42|l 
and (|43|l is rather simple: if we want to get a self- 
consistent distribution which coincides with a solution 
of the S-^ equation for a given zero perveance poten- 
tial, the simplest way it is to calculate the space charge 
potential for this frozen distribution through the Poisson 
equation, and then compensate the external potential ex- 
actly for that. This is what we did in our example where 
the gaussian solution is the fundamental state of a har- 
monic oscillator: we finally got a total potential which is 
vo{s) = v{s) + fe(s) — (namely that of a simple har- 
monic oscillator), and an energy value which coincides 
with the first eigenvalue. In other words, if you want a 
gaussian transverse distribution you should not simply 
turn on a bare harmonic potential s^: you should rather 
teleologically compensate for the space charge by using 
the potential Ue(s). 



B. Student transverse distributions 

If the halo consists in the fact that large deviations 
from the beam axis are possible, a new idea is to suppose 
that the the stationary transverse distribution is different 
from the gaussian distribution (|51|) introduced in the Sec- 
tion [IV^ To this end we will introduce in the following 
a family of distributions which decay with the distance 
from the axis only with a power law. 

Let us consider the following family of univariate, two- 
parameters probability laws S(i^, a^) characterized by the 
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0.^ 




FIG. 2: The Gauss pdf Af{0, 1) (dashed line) compared with 
the E(2,2) (thick hne) and the E(10, 12) (thin line). The 
flexes of the three curves coincide. Apparently the tails of the 
Student laws are much longer. 



following pdf 's 



'£+1 ' 



(55) 



which apparently are symmetric functions with the mode 
in a; = and two flexes inx = zta/^/i/ + 2. All these laws 
are centered at the median. In particular a plays just 
the role of a scale parameter, while v rules the power 
decay of the tails: for large x the tails go as with 
v + 1 > 1. For a comparison with a Gauss law Af{0, cr'^) 
see FIG. 121 Remark that when v grows larger and larger, 
the difference between the two pdf 's becomes smaller and 
smaller. It is typical of the laws £(1/, a^) that they have 
(finite) momenta of order k only if the condition k < v 
is verified; hence for v <2 there is no variance, while for 
v < \ not even the expectation is defined. On the other 
hand when v > 2 the variance of S(j^, a?') exists and is 



(56) 



It will be useful to remark that the laws E(l, a'^) are the 
well-known Cauchy laws C(a) with pdf 

fix) = -- 



while the laws S(n, n) with n : 
i-Student laws S{n) with pdf 



1,2,... are the classical 



We will then refer to a?) as generalized Student laws 
since they are just Student laws with a continuous param- 
eter V > Q and a scale parameter a. For v > 2 variances 
exist and we are then entitled to standardize our laws: 
indeed from every E(z^, {v—2)(j'^) with a? = {v—2)a'^ 



has variance ct^, and the standard (with unit variance) 
generalized Student laws are S(i/, v — 2). 

In order to describe the beam we will also introduce the 
bivariate, circularly symmetric Student laws Y^2{i'tO?) 
with pdf 



.f{x,y) 



27r 



(x2 



(57) 



Its marginal laws are both I](i^, a^) and non-correlated, 
albeit not independent (as in the case of the circularly 
symmetric gaussian bivariate laws) . The total beam dis- 
tribution will then be 



H 



z\ (58) 



where H{z) is the Heaviside function. In the descrip- 
tion of a beam in an accelerator it is realistic to suppose 
that the transverse distribution is endowed with a finite 
variance. Hence we will look for distributions (|58|) with 
v > 2. On the other hand this will correspond to suppose 
that in our model the transverse Student laws should not 
be radically different from a Gaussian: in fact the halo 
is in some sense an effect which is small when compared 
with the total beam. From this standpoint the family 
of laws I](z^, a^) has also the advantage that we can fine 
tune the parameters v, a in order to get the right dis- 
tance from the gaussian laws (this would not be possible 
if we adopted stable laws; see subsequent Section IVAp . 
With this hypothesis in mind we will limit our present 
considerations to the case z/ > 2 so that the transverse 
marginals of H58|l will have a finite variance a^. Then 
from (|56|l we choose = (ly — 2)0^ and write (|58|l as 



[{v-2)a^]^ 



H[k-\z\ 



2ttL [2^2 _^y2 ^ (^j^ _ 2)cr2] -T- \2 



(59) 



Passing to cylindrical random variables we then have 

^2 



p(r, ip,z) =r 

namely 

p{r,ip,z) 



[(^-2)a2]^ ^/L 
2t:L [r2 + (i.-2)a2]^ ^2 



1 



21^ 



aV2 crV2 v~2 



(i.- 2)a2 



2ttL 



so that finally with the shorthand notation 

s^/2 



the dimensionless, normalized radial distribution is 



2v 



1 



2(l + z2)^ 



(60) 



(61) 
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Here we adopt the dimensional constants 
„ 2 

A = cr\/2 



V 



(62) 



where is the variance of our Student laws. We can 
now use the relations (|^ . and in order to get 
the potentials which have (|59ll as stationary distribution: 
first of all the space charge potential produced by H59|) 
has the form 



v{s) 



2z- 



2F1 



1 



+ log + C + V 



(63) 



where 2Fi{a,b;c;w) is a hypergeometric function and 
tjj{w) — T'{w)/r{w) is the logarithmic derivative of the 
Euler Gamma function (digamma function). On the 
other hand, by choosing /3 = 2 + to put the po- 
tential energies to zero in the origin, we get the control 
potential for zero perveance 



wo(s) = 



iy + 2 z2(4z2 + iy+10) 

1^-2 2(1 + ^2)2 



(64) 



and hence the external potential required to keep a trans- 
verse student distribution 'S2{v, (v ~ 2)a'^) with a given 



variance a is 



Ve{s) 



2(l + z2)2 

■2z-^ 
2-ri 



-logz^ + C + V - 



(65) 



Formulas H63() . H64|) and (|65(l give the self-consistent po- 
tentials associated with the beam distribution H59(l which 
is transversally a Student S2(z^, {i^ — 2)a^). In the FIG. El 
we can see an example of the control potential vo{s) for 
a particular value of the parameter v, together with its 
limit behaviors 



vo{s) 
vo{s) 



{iy + 2){i^- 



M,2, 



(z.-2)2 

(^ + 2)^ , , 



4s2 



(s ^ 0+) (66) 
(s^+oo) (67) 



Now this results must be compared with the similar re- 
sults (|52() . H53() and H54|l associated to a transversally 
gaussian distribution. We will choose the gaussian pa- 
rameters in such a way that the behavior near the beam 
axis be similar to (|65|l . namely (with (3 = 27) 



w(s) 



276" 



-27SV2 



7 



{iy + 2){iy+10) 
{V-2Y ■ 



First of all in the FIG. 01 we compare the space charge 
potential produced by both a Student and Gauss trans- 
verse distribution: remark as for the chosen parameter 
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FIG. 3: The control potential vo{s) 16411 for a Student trans- 
verse distribution E2(22, 20o-^). Also displayed are the value 
of /3 = 2.4 (the limit value of vq for large s, thin line) and the 
behaviors for small and large s H66|l and l|67|l (dashed lines). 



values {v — 22, ^ — 20) the two potentials look particu- 
larly similar. In fact, given the asymptotic behavior of 
the hypergeometric function in H63|) and of the exponen- 
tial integral in (|52|l . for s +00 both potentials behave 
as — ^logs. On the other hand we immediately see from 
FIG.[Slthat the control potentials for zero perveance vq{s) 
behave differently when we move away from the beam 
axis; beyond a distance of about r ~ 2ct the two curves 
are different: while in the Gaussian case the potential 
diverges as s'^, in the Student case it goes to the con- 
stant value P as quickly as s~2. Of course this difference 
fades away when v grows larger and larger; that points 
to the fact that the principal difference between the two 
cases can be confined in a region that can be made as far 
removed from the beam core as we want by a suitable 
choice of i'. Finally in the FIG. Owe compare the total 
external potentials needed to keep the transverse beam 
respectively in a Student and in a Gauss distribution. We 
then see that for large s (far away from the beam core, 
while in the Gauss case the total external potential grows 
with s as s2 + ^ log s, in the Student case this potential 
only grows as ^ log s. In any case, even if the potential 
near the beam axis is harmonic, deviations from this be- 
havior in a region removed form the core can produce a 
deformation of the distribution from the gaussian to the 
Student. 



C. Estimating the emittance 

If u{r) is a self-consistent, cylindrically symmetric so- 
lution of H36|l and (|37|l the position probability density 
p(r, If, z) in cylindrical coordinates will have the form 



p{r,ip,z) = 



1 



2(r) 0<ip<2TT, 



2 ^ ^ 2 ' 



27ri 



otherwise. 



In order to estimate the emittance we need to calculate 
mean values of positions x and momenta Px along one 
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FIG. 4: The space charge potentials v{s) ib'Sti and 1521 re- 
spectively for a Student (solid line) transverse distribution 
E2(22, 20(T^), and for a Gauss (dashed line) distribution. The 
dimensionless perveance here is ^ = 20. 




FIG. 5: The control potential vo{s) H64|l of a Student 
E2(22,20cr^) (solid line; see FIG. is here compared with 
that of a Gauss distribution (dashed line) which shows the 
same behavior near the beam axis. 




FIG. 6; The total external potential Ve{s) 16511 that should 
be applied to get a stationary Student transverse distribution 
E2(22, 20(T^) (solid line), compared with that l|54^ needed for 
a Gauss distribution (dashed line). 



transverse direction, but we should remember that in SM 
we have neither a distribution in the phase space, nor an 
operator formahsm. The momentum and its distribu- 
tion should then be recovered from the velocity fields Q 
and Q where ~ since we are dealing with stationary 
states with v = - only the osmotic part is non zero 
so that 

p = a— = 2a \' . 

p u[r) 

By supposing now to choose the v of our Student laws so 
that the following integrals exist, we then have 







dz 



dz 

T 



27T 



COS Lp 

2tt 



dip 



+ 00 



r^M^(r) dr = 







cos^ p 

271 



dp 



+ 00 



r u (r) dr 



(Px) 

{pD 



1/ .3.,. 

2 70 

2a — 

Jo L jQ 



r' u (r) dr 
L 



cos </? 



27r 



dp 



ru{r)u'{r) dr = 



cos^ p 
2n 



^ f dz 
Aa^ - 
Jq ^ JO 

r+oc 

2c? \ ru''^{r)dr 
Jo 



dip 



ru' (r) dr 



so that the standard deviations (uncertainties) are 



Ax = 




r^u'^{r) dr 



(68) 



Apx 



-\-oo 



'2a2 I ru'^{r)dr (69) 



and the position-momentum covariance is 

C = {xpx) - {x){px) = (xpx) 

^ I dz I con p I 2 t / \ / \ 1 

^ 2a — — dp / r u [r)u[r) dr 

Jo L 7o 27r Jq 

p + QC 

— a I r'^u' {r)u{r) dr 
Jo 

In a previous paper |3] we adopted the uncertainty prod- 
uct Ace • Apx as a measure of the r.m.s. emittance. As an 
example let us suppose again that our wave function has 
the form Uo(r) for the harmonic oscillator without space 
charge given in (|45(l . We then have 



Aa; • Apx 



^C^a. 



(70) 



This allows two remarks: first, a plays also the role of 
a measure of the emittance and hence - as suggested in 
a previous paper - its value must be linked to the 
number of particles in the beam; second, the position- 
momentum correlation coefficient of a Gaussian beam is 

C 



Ax ■ Apx 



-1 
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as it was predictable, since in SM the relation between 
position and momentum for the wave function H45|) is 
linear and negative. 

In other models the transverse r.m.s. emittance is cal- 
culated by means of the quantity ^ IS.x^ Ap^ — . In 
the KV distribution, since momentum and position are 
uncorrelated and (x) = (p^:) = 0, this estimate be- 
comes yj{x'^){p^) — (xpj.)^. In the SM model, on the con- 
trary, this is not a good choice: in fact we have shown, 
at least in our simple example, that x and are far 
to be uncorrelated, and that as a consequence of H70() 
y/ Ax^ Ap^ — (72 becomes exactly zero. Apparently it is 
not realistic to take this value as a good estimate of the 
emittance. On the other hand, for the same gaussian 
example, the value of the uncertainty product Ax Ap^ 
is just a which we assume to be a good candidate for 
the value of the emittance. On the other hand it is easy 
to calculate the same uncertainty product for a Student 
distribution E2(^, (z^ — 2)cr^) with dimensionless radial 
distribution Ht)l|l and variance a^: in fact a straight ap- 
plication of (|68|l and H69() brings to the following result 



Aa; • Apx — a 



(71) 



Of course, as it is already clear, this value converges to 
the Gaussian case for large while becomes larger and 
larger for small v values when the shape of the distribu- 
tion moves away from the Gaussian case. 



D. Weighing the tails 

We can finally compare the length of the tails of Gauss 
and Student distribution in order do assess the possible 
halo formation in the second case. Let us consider the 
probability 



P(c) 



+ 00 



(72) 



of being beyond a distance ccr (cr^ being the variance) 
away from the beam axis, and calculate this quantity in 
our two cases. From the Gaussian distribution we have 
from ESI that 



practically no particle beyond lOcr in the Gaussian case, 
but about 10'^ particle per meter for & v — 22 Student 
distribution, and as much as 10^ for a.v — IQ value. It is 
worthwhile to remember at this point that we got about 
the same number of particles gone astray in our self- 
consistent numerical solutions for a dimensionless per- 
veance of about ^ = 20 in one of our previous paper Jll . 



V. LEVY-STUDENT PROCESSES 

In our context the Student laws S(i^, a?) are important 
not only because they promise to better describe the halo 
by means of their longer tails with respect to usual Gaus- 
sian distributions; in fact they constitute an important 
family of Levy infinitely divisible (i.d.) laws. At present 
there is a lot of interest about non-Gaussian Levy laws 
in several fields of research (see for example ^3 and 
references quoted therein), but this interest is mostly con- 
fined to the stable laws which are in fact an important 
sub-family of the i.d. laws. The fundamental character 
of the i.d. laws can be better understood from two differ- 
ent, but strictly correlated standpoints: on the one hand 
the i.d. laws constitute the more general form of possi- 
ble limit laws for the generalized Central Limit Theorem; 
on the other they constitute the class of all the laws of 
the increments for every stationary, stochastically con- 
tinuous, independent increments process (Levy process). 
These important results (which are briefly discussed in 
Appendix 1X1 and Appendix IbI) have been achieved by P. 
Levy, A.Ya. Khintchin, A. Kolmogorov and other mathe- 
maticians from the mid 30's to the mid 40's of the XXth 
century, but their relevance for the applications has been 
recognized only in more recent years. One of the charac- 
teristics of a non-Gaussian Levy process is to have tra- 
jectories with moving discontinuities (think to the trajec- 
tories of a typical Poisson process contrasted with those 
of a Gaussian Wiener process), and we propose here to 
describe the trajectories of the particle beam by means of 
a Levy-Student process whose discontinuities can possi- 
bly account for the relatively rare escape of particles from 
the beam core. For the sake of simplicity we will limit 
ourselves in the following to the case of 1-DIM systems 
representing one single transverse coordinate of our par- 
ticle beam. 



P{c) = e-^'/^ (73) 
while in the Student case from (|59|l we get 



P{c) 



v-2 



(74) 



Now for c = 10 the Gaussian value is about 1.9 x 10^^^, 
while with = 10 the Student value is about 2.2 x 10~®, 
and with 1/ = 22 the value is 2.8 x 10~^. This means 
that for M = 10^-^ particle per meter of beam, we find 



A. The Student i.d. laws 

The ch.f.'s of the laws E(j^, a^), namely the Fourier 
transform of the densities (|55|l . are 



|qk|^K^(|qk|) 



(75) 



where Ka(z) is a modified Bessel function. The typical 
form of these ch.f's (contrasted with the Gauss ch.f.) is 
shown in FIG. [7| Remark that, since a; is a length, the 
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FIG. 7: Typical ch.f. of a Student law E(2, 2) (solid line) 
compared with the ch.f. of a standard Gauss law N{Q, 1) 
(dashed line). 



only for i — s = T, while for t — s ^ T it is no more 
the ch.f. of a Y,{v,a'^) law. Hence we explicitly know 
the form of the increment law only at the time scale T: 
we know the ch.f. - namely everything we theoretically 
need - at every time, but we do not have the explicit 
inverse Fourier transform, and we also know that the 
laws are no more in the family S(i^, a^). This problem is 
tempered by the remark that the situation is not better 
for general stable laws: even in this case, in fact, we do 
not know the explicit forms of the increment laws not 
even for one time scale (they are known only in a few 
precious instances). The unique advantage in the stable 
case being the fact that all along the time evolution the 
increment laws remain of the same type, which is not the 
case for i.d. non stable laws. 



B. The Levy— Student process 



ch.f. variable k has the dimensions of wave number (in- 
verse of a length). These laws are i.d. but in general are 
not stable, the unique stable laws among them being the 
Cauchy laws C{a) = E(l,a^). For v > 2 the Student 
distributions belong to the domain of attraction of the 
Gauss law since they have a finite variance. For v < 2 
the variance diverges, but it is possible to prove that this 
notwithstanding they still belong to the Gaussian domain 
also for V = 2: albeit this derives from a well known gen- 
eral result a simple proof for our particular case will 
be given in a subsequent paper. On the other hand for 
u < 2 the Student laws are attracted by non-Gaussian 
stable laws characterized by the same value of the pa- 
rameter V. 

The fact that the Student laws are i.d. - which in itself 
is not at all a trivial result proved in steps only in the 
70's and 80's 129| - shows two kinds of advantages with 
respect to more common stable laws: 

• no stable, non-Gaussian law can have a finite vari- 
ance, while all Student laws with v > 2 do have 
a finite variance; this is important since it is not 
realistic to suppose that empirical distributions (in 
particular for the particle beams) have infinite vari- 
ances, but that notwithstanding we will not be 
obliged to resort to handmade modification (for in- 
stance truncated Levy distributions) as in the case 
of stable distributions |lflj|; 

• the asymptotic behavior of stable, non-Gaussian 
laws is proportional to |a;|-"--^ with a < 2 [l8l |. 
while the asymptotic behavior of the Student laws 
is |a;|-'^-"'^ with v > Q] this allows the Student 
laws - but not the stable laws - to continuously 
go throughout all the gamut of decay speeds to ap- 
proximate in a fine tuning the Gaussian behavior 
as well as we want. 

The principal drawback for not being stable is in the 
subsequent definition of the Levy-Student process. In 
fact the ch.f. of the process <p(k, t) coincides with H75|) 



A Levy process defined by the ch.f. (|75|l will be called 
in the following a Levy-Student process. Taking into ac- 
count (jB3|l and (|75|l we can now state that the transition 
pdf of a Levy-Student process is 



p{x,t\y,s) 



1 

2^ 



+00 



\aK\2KiL{\aK\) 



2^r(|) 



dK (76) 



where the improper integral is always convergent since 
the asymptotic behavior of the ch.f. is 



(P{k) 



12-K 



[1 + 0(1^^1-1)] 



2§r(|) 



-CX) 



In principle (|76() should be enough to calculate every- 
thing of our process, but in practice this is an integral 
that must be treated numerically, but for a few particu- 
lar cases that will be discussed in a subsequent paper. On 
the other hand even to produce simulation for the tra- 
jectories of our process we should have some simple ex- 
pression for the transition pdf. At least for this last task, 
however, we can exploit the fact that when t—s = T the 
expression (|76|l can be exactly calculated and coincides 
with the pdf (|55|) of a Student I](j/, a^) (remark that even 
this is not possible for the typical non-Gaussian stable 
process). This means that we can produce sample tra- 
jectories by taking T as the fundamental step of our nu- 
merical simulation. In other words we will simulate the 
sample paths of a process whose increments are exactly 
Student distributed when observed at the (otherwise ar- 
bitrary) time scale T. To give a look to these trajectories 
we produced a simplified model which simulates the so- 
lutions of the following two SDK's 

dX{t) = v{X{t),t)dt + dW{t) (77) 
dY{t) = v(Y{t),t)dt + dS{t) (78) 

where W{t) is a Wiener process, while S{t) is a Levy- 
Student process. We also fixed the velocity field v{x,t) 
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FIG. 8: The pdf's of the increments for the Gaussian pro- 
cesses (SDE (1771 . dashed line; a ~ 0.53) and for the Levy- 
Student process with law E(4, 1) (SDE 17811 . solid line; a ~ 
0.71). The parameters are chosen so that the two pdf's have 
the same modal values and similar shapes. 
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FIG. 9: Typical trajectory of a stationary, Gaussian 
(Ornstein-Uhlenbeck) process X(t) (see SDE l|Z3). To com- 
pare it with the Student trajectory, the vertical scale has been 
set equal to that of FIG.ITni 

in a suitable way: it will not depend on time t, and its 
value is (for given 6 > and q > 0) 

v{x) — ~bxH{q~ \x\) 

where H is the Heaviside function. This flux will at- 
tract the trajectory toward the origin when |a;| < 
and will allow the movement to be completely free for 
\x\ > q. The forms of the typical pdf's used in our sim- 
ulations are shown in FIG. |H| In a simplified model for 
a collimated beam this will then produce a stationary, 
Ornstein-Uhlenbeck process for the SDE H77() if the in- 
tensity of the Gaussian noise is not too large. The pro- 
cess solution of the SDE H78|l will instead have different 
characteristics. Let us suppose to fix the ideas that the 
two parameters defining the velocity field are b — 0.35 
and q = 10. The FIG. El displays a typical trajectory 
of a 10^ steps solution X{t) of H77|) when the variance 
of the Gaussian distributed increments is = 0.28. In 



FIG. 10: Typical trajectory of a stationary, Student process 
Y{t) (see EDS (EHJ, v = 4 and a = 1). 



our simplified 1-DIM model of the transverse dynamics 
of a particle beam this means that the trajectories al- 
ways stay inside the beam core. Let us then take as law 
for the increments of (|78|) a Student distribution I](4, 1): 
its pdf looks not very different from that of the previ- 
ous Gaussian distribution, as the FIG. |H1 clearly show. 
That notwithstanding the process Y{t) differs in several 
respects from X{t). Indeed not only the typical trajec- 
tory displayed in FIG. llOl shows a wider dispersion of its 
values and a few larger spikes. The principal difference is 
rather in the fact that while the trajectories of X{t) show 
a remarkable stability in their statistical behavior, the 
paths of Y{t) have the propensity to make occasional ex- 
cursions far away from the beam core (see FIG. Ill|l. and 
seldom they also definitely drift away from the core (see 
FIG. I12|I . This depends of course on the mentioned prop- 
erties of the trajectories of a non-Gaussian Levy process, 
and in particular on the fact that they are only stochas- 
tically, and not pathwise continuous, namely that they 
contain occasional jumps. The frequency and the size of 
these jumps can also be fine tuned by suitably choosing 
the values of the parameters of the law E(i/, a^) of the 
increments. It is this feature of a Levy-Student process 
that suggests to adopt this model to describe the rare 
escape of particles away from the beam core. 



VI. CONCLUSIONS 

In the previous sections we have introduced the Stu- 
dent laws in our SM model for the particle beam dy- 
namics first of all in order to make use of their features 
depending on their enhanced variance. In particular we 
have shown that the longer tails with respect to the sim- 
ilar Gaussian distributions can help to account for the 
finding of a larger than expected number of particles re- 
moved far away from the beam core. 

It should be remarked, however, that all along the Sec- 
tion ^3 our processes were Gaussian processes since the 
underlying SDE Q is still powered by a Brownian noise. 
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FIG. 11: Occasional trajectory of a stationary, Student pro- 
cess Y{t) (see EDS 17811 . f = 4) with a temporary excursion 
out of the core. 




FIG. 12; Rare, but possible trajectory of a stationary. Stu- 
dent process Y{t) (see EDS 17811 . i/ = 4): here the particle 
definitely drifts away from the core. 



for a process produced by a Levy-Student noise (a pro- 
cess only stochastically continuous: trajectories can have 
jumps) occasionally the jump is large enough to put the 
particle out of the stream. Of course the frequency and 
the size of these jumps depend on the parameters v and 
a of the process: the jumps tend to be smaller and less 
frequent when the S(i^, a^) distributions approximate a 
Gaussian law. In our opinion it would be very interest- 
ing to explore the possibility that the processes under- 
lying the intense beam dynamics be ruled by some sort 
of Levy-Student noise rather than by the usual Gaus- 
sian noise. It is then important to point out that a few 
numerical evidences ,20j begin to emerge which confirm 
this conjecture. 

These remarks point to several research directions. 
First of all it is important to better study the Levy- 
Student process in itself: for example a knowledge of the 
Levy-Khintchin functions of the Student laws would be 
relevant to the fine tuning of the frequency and the size 
of the trajectory jumps. On the other hand even the dif- 
ferential form of its Chapman-Kolmogorov equation 's^l 
would be instrumental to discuss the time evolution of 
the process. Then it must be remarked that at present 
we have just defined the Levy-Student process, but we 
added no dynamics: it is as if we have the Wiener pro- 
cess, but no Stochastic Mechanics or any other dynamical 
model added to this kinematics. In other words we need 
to build a new generalized SM for the Levy-Student pro- 
cesses. Finally it would be important at this point to 
have empirical or numerical data able to corroborate the 
hypothesis that the increments of the transverse variables 
of a beam are in fact distributed according to a Student 
law, rather than according to the usual Gaussian law. 
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This is true even in the Section flV BI where we first intro- 
duced the Student laws IjSSI) as stationary distributions 
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APPENDIX A: INFINITELY DIVISIBLE AND 
STABLE LAWS 

The relevant mathematical concepts used in this paper 
are better discussed in the framework of the theory of the 
addition of independent random variables (r.v.): for more 
details see |^ [S^, ■ In the following we will describe 
the law C oi a, r.v. X by giving her characteristic function 
(ch.f.) 

^{k) = E(e-^) 

where E(-) is the expectation under the law £. When £ 
has a pdf f{x), then (p^n) is just its Fourier transform. It 
is well known that the law C of the sum of n independent 
r.v.'s with laws £i, . . . , Cn has a ch.f. which is the product 
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of the ch.f.'s of the component laws: 

IP{k) ^ IPi{k) ■ . . . ■ ipn{K) (Al) 

On the other hand we say that a law C is decomposed 
in the laws Ci, . . . ,Cn when its ch.f. can be written as 
a product IjAljl of the ch.f.'s of its components. This 
already allows us to introduce two fundamental concepts: 
a law C with ch.f. ip is said to be i.d. when for every n 
there is a law Cn with ch.f. ipn such that ip — ip"^. In other 
words this means that for every n a r.v. X with law C 
can always be decomposed in the sum of n independent 
r.v.'s all with the same law (identically distributed). 
Remark, however, that in general the laws £„ are not 
of the same type as C. Let us remember here that we 
say that two laws are of the same type when we get one 
from the other by means of a centering and a rescaling; 
in other words, if in a ch.f., then all the ch.f.'s of the 
same type have the form e^°''^(p{bK) for every a and b > 0. 
For instance all the Gaussian laws A/'(/i, cr^) belong to the 
same (Gaussian) type; on the contrary the Poisson laws 
7'(A) with different values of A do not belong to the same 
type. Now, a law £ is said to be stable when it is i.d. 
and the component laws are of the same type as £. More 
precisely a ch.f. (^(k) is stable when for every b,b' > 
there exist a and c such that 

^(ck) = e""'ip{bK)ip{b'K) . 

As an example: the Gaussian and the Cauchy laws are 
stable; the Poisson laws are instead only i.d. The fami- 
lies of i.d. and stable laws are completely characterized: 
in fact the celebrated Levy-Khintchin formula gives the 
more general form for the ch.f.'s of these two classes; 
however, while in the case of the stable laws these ch.f.'s 
(albeit not in general the laws themselves) are explicitly 
known in terms of elementary functions, for the i.d. laws 
the ch.f.'s are given through an integral containing a func- 
tion L(x) {Levy function) associated to every particular 
law. But for a few classical cases the Levy functions of 
the i.d. laws are not known. 



APPENDIX B: CENTRAL LIMIT THEOREM 
AND LEVY PROCESSES 

Let us consider the sequence of r.v.'s Xn^k with n e N 
and k = 1, . . . , n with A„_i, . . . , A„^„ independent for 
every n. The modern formulation of the Central Limit 
Problem asks to find the more general laws which are 
limits of the laws of the consecutive sums 

n 

Sn = ^ Xn,k (Bl) 
fc=l 

Remark that these sums generalize the usual partial sums 
of the classical Central Limit Theorem in that: when we 
go from Sn to, say, Sn+i, the first n terms do not in gen- 
eral remain the same: for example Xn.i does not coincide 



with Xn+i.i- Under very general technical conditions the 
Central Limit Theorem now states that the family of all 
the limit laws of the consecutive sums (jBl|) coincides with 
the family of i.d. laws. The stable laws come into play 
only when we specialize the form of our consecutive sums: 
when we have 

_ Xk bn 

-^n,k — 

an n 

where a„ and 6„ are sequences of numbers, and Xk are 
independent r.v.'s, the consecutive sums take the form of 
the usual normed sums (centered and rescaled sums of 
independent r.v.'s) 

c* " 
Sn^^-bn, 5:=^Afc. (B2) 

k—1 

Then, if the Xk are also identically distributed, the family 
of the limit laws of the normed sums ljB2p coincides with 
the family of the stable laws. The classical (Gaussian) 
Central Limit Theorem is an example of convergence to- 
ward a stable law; on the other hand the Poisson Theo- 
rem (convergence of Binomial laws toward Poisson laws) 
is an example of convergence toward an i.d. law. Every 
stable law has its own domain of attraction, namely the 
set of laws attracted by it in the sense of the convergence 
of normed sums ljB2|) of independent r.v.'s all distributed 
as the attracted law. It can be proved that all the laws 
with finite variance are in the domain of attraction of the 
Gauss law, and that a law can be attracted by a non- 
Gaussian stable law only if it has infinite variance. 

The general formulation of the Central Limit Theo- 
rem is strictly connected to the definition of the pro- 
cesses with independent increments (decomposable pro- 
cesses). It is apparent in fact that if the increments 
AX{t) = X{t + At) — X{t) for non superposed inter- 
vals are independent, the previous forms of the Central 
Limit Theorem imply that the laws of the increments 
must be i.d. laws. Moreover, since the decomposable 
process are also Markov processes, the laws of the in- 
crements are also all that is needed to completely define 
them. If a decomposable processes X{t) is stationary 
(namely the law oi X{t + s) — X(s) does not depend on 
s) and stochastically continuous (namely for every t we 
have X{t + At) - A(i) ^ in probabihty when At 0) 
we will call it a Levy process. Remark that a Poisson pro- 
cess is a Levy process since, despite its discontinuities, it 
is stochastically continuous. In fact these discontinuities 
do not impair the stochastic continuity of the process be- 
cause they are moving (as opposed to fixed) discontinu- 
ities. On the other hand it is possible to prove that only 
the Gaussian Levy processes (for example the Wiener, or 
the Ornstein-Uhlenbeck processes) are pathwise contin- 
uous, namely: almost every sample path is everywhere 
continuous (there are not even moving discontinuities). 
Now, if (p{k) is the ch.f. of an i.d. law and T is a suit- 
able time constant, it is possible to prove that [(^(k)]^'/-^ 
is the ch.f. of the increments AX{t) of a Levy process. 
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Hence, if the process has a pdf, the stationary transition 
pdf is 

pix, t\y, s) = — PV / e-(-^) [ip{K)] ^ dK (B3) 

SO that, at least in principle, we know all that is needed 
to define the process. 

The sample paths of a Levy process are also well char- 
acterized: it is possible in fact to prove that almost all 
trajectories are bounded and are continuous with the ex- 
ception of a countable set of moving jumps (first kind 



discontinuities). Then, let us suppose that Lt{x) is the 
Levy-Khintchin function of the i.d. law of the increment 
X{s + t) — X{s): if iyt{x) is the random number of the 
jumps in [s, s + t) of height in absolute value larger than 
x > 0, it is possible to prove that 

\Ltix)\=Ei,yt{x)) 

so that the Levy-Khintchin function of an i.d. law plays 
also the role of a measure of the frequency and height of 
the trajectory jumps. 
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